Freezing of Triangulations 



J. -P. Kownacki 



Laboratoire de Physique Theorique et Modelisation, 
CNRS-Universite de Cergy-Pontoise, UMR 8089 
5, mail Gay-Lussac, Neuville-sur-Oise, 95031 Cergy-Pontoise Cedex, France 
E-mail : kownacki@ptm.u-cergy.fr 



Abstract 

Zero temperature dynamics of two dimensional triangulations of a torus with curvature energy is described. 
Numerical simulations strongly suggest that the model get frozen in mctastable states, made of topological 
defects on flat surfaces, that group into clusters of same topological charge. It is conjectured that freezing is 
related to high temperature structure of baby universes. 
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1 Introduction 



Two dimensional dynamical triangulations models have been largely studied in the context of euclidean 
quantum gravity over the past decades (for reviews, see [1]). These models allow lattice regularization 
where the cut-off is the length of the links. In view of 2d quantum gravity formulation, interesting 
features are essentially thermodynamical limits and critical static properties, as they are key features 
for continuum limit of lattice formulation. Great analytical progresses have been made but numerical 
simulations are particulary suited to these models and allow further investigations in remaining open 
problems. Square curvature terms in the action/hamiltonian is such a problem that still lacks an 
analytical solution. 

However, as random surfaces models, dynamical triangulations can be used in several other domains 
like real membranes, foams, 2d-liquid, defects on surfaces, etc. [2] [3]. Dynamical triangulations models 
have the nice feature to be very simple to formulate - purely topological models - and to display very 
interesting physical behaviours. 

Recently, Sherrington et al. have studied the dynamics of such a model with curvature energy. 
They discovered behaviours characteristic of strong glass formers [4]. They underlined the role of 
topological constraints in this model without quench disorder, and explained glassy behaviour as 
dynamics of topological defects at low temperature. Their aim was to study a simple model of 
supercooled liquid, but their results are not limited to it. In particular, dynamical properties are 
essential in numerical simulations, even if static properties are studied - as in 2D quantum gravity - as 
Monte-Carlo techniques consist in performing random walks in configuration space. In such a context, 
a glassy behaviour with exponentially long equilibration times can mask true equilibrium properties. 

In this paper, the work described above is extended to zero temperature dynamics of a square 
curvature energy model. In fact, zero temperature limit, where no thermally activated process is 
allowed, is crucial to understand, as it gives more informations on the energy landscape and, thereby, 
on the whole low temperature dynamics. 

2 The model - Static properties 
2.1 Definition 

The set of all regular two dimensional triangulations [5] with N vertices, toroidal topology and no 
boundaries is considered. Regular means triangulations without tadpoles - vertex connected to itself 
- and self-energy diagrams - two vertices connected by two links. Moreover, each vertex is connected 
to at least three neighbours. This set is noted T^ orus . Topological constraints imply that the number 
of vertices N, the number of links N and the number of triangles N t are not independent. The first 
constraint sets the Euler characteristic x to zero for a torus, x is a topological invariant depending 
on the number of classes of nontrivial loops on the surface - two for a torus. For a triangulation, it 
reads [6]: 

X = N-N l + N t 

so that, for a torus, 

N - N + N t = (1) 
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The second constraint is local: each triangle is made of three links and, if there are no boundaries, 
each link belongs to two triangles. So, 

2N t = 3N t (2) 

Combining relations (1) and (2) gives 

N t = 3N , N t = 2N 

In the model described here, each triangulation T £ 2 1 * orus i s endowed with one - and only one - 
geometry. This is achieved by setting the length of all links to unity and imposing on metric to be flat 
inside each triangle of T. Such geometries are called piecewise linear manifolds. It should be noticed 
that triangulations are not embedded in any larger space, so that geometries are internal. 

A random surfaces model is built from this ensemble of surfaces with fixed topology but different 
geometries. The method consists to assign to each triangulation/geometry an energy E(T) which 
governs dynamics and equilibrium. At equilibrium, each triangulation T is weighted by a Boltzman 
exponential exp (—(3E(T)) with f5 the inverse temperature. In the model, E(T) depends on Gaussian 
- internal - curvature. As geometry inside each triangle is flat, curvature is concentrated at vertices. 
Following Regge calculus [7], an elementary surface dsi = qi is defined at each vertex (i) connected to 
qi neighbours (see figure 1). Such a vertex is called a ft-vertex. 



Figure 1: A triangulation. 
Local curvature at vertex (i) is qiRi, with 

r-> (6-gi) 
Ui — 

Qi 

Ri measures local deviation from flatness (qi = 6 and R4 = 0). However, E{T) cannot be the total 
curvature ^ QiRi because it is a topological invariant, as stated by the Gauss-Bonnet theorem, and 
would be a constant in the model. The choice for the energy is 

E(T) = ]T ( ft - 6) 2 (3) 

i 

It is not strictly speaking equal to the total square curvature 
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but (3) and (4) are expected to give similar results, at least for small curvature, because 
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In both cases, surfaces are flattened by curvature energy at low temperature. However, expression (3) 
is more symetric with respect to q% = 6. Moreover, for surfaces containing only 5-, 6- and 7-vertices, 
the energy can be rewritten 



where is the number of /c-vertices. As will be seen later, these surfaces play a crucial role at low 
temperature. 

Finally, the partition function of the model reads 



where C(T) is a symetry factor needed to avoid overcounting symetric triangulations. For large N, 
C(T) = 1 for almost all triangulations. 

2.2 Equilibrium Properties 

At infinite temperature (/? = 0), the model has been solved [8] [9]. Typical surfaces are characterized 
by hierarchical structures of baby universes - bubbles growing on surfaces and linked by very small 
necks. The distribution of their size is known [10]. As a consequence, typical surfaces are fractal 
with Haussdorf dimension dn = 4 [11]. The distribution of vertices with n neighbours P(n) has been 
calculated [12] (see figure 2 ) and the variance is fi2 = 10.5, meaning a wide range of vertices on the 
surfaces. Note that the model follows a modified Aboav's law [13]. 
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Figure 2: Equilibrium distribution P(n) at f3 = 0. 
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At finite temperature, the model has not yet been solved and static properties have been studied 
numerically by many authors [14] [15]. However, a very similar model has been solved analytically, 
using matrix model techniques [16]. 

When temperature is lowered, curvature energy favours flattening of surfaces. The ground state of 
the model is the hexagonal lattice with % = 6 at each vertex (figure 3). It should be emphasized that 



Figure 3: The hexagonal lattice, 
the ground state does not depend on the precise form of the energy. In fact, E(T) can be rewritten 

E{ 7) = E (V + 36-12^) 

i 

= E (V) + 36N - 2AN i 



(V) - 36iV 



so that, at given TV, E(T) is equivalent to J2i(Qi ~ Qo) 2 , whatever the value of q Q . On the other hand, 
the mean neighbours number 

<«> = n E* 

i 

= T7 2 N l 
N 

= 6 

is completely determined by topological constraints (1) and (2). So, the ground state is the configu- 
ration that minimizes J2i Qi with < % >= 6. 

Here are summarized results obtained by numerical simulations, done as preliminary work in this 
paper, for sizes up to iV = 65025. Starting from hexagonal lattices, the system is initially equilibrated 
with typically 1000 sweeps, and simulations are carried with about 3000 sweeps. For details about the 
rules of simulations, see paragraph "dynamics" below. As can be seen in figure 4, the distribution P(n) 
is more and more peaked at n = 6 as temperature is lowered, meaning that surfaces contain essentially 
6-vertices. Then, 5- and 7-vertices can be seen as topological defects in flat surfaces, whose proportion 
vanishes as temperature is lowered. It should be noticed that, when only 5- , 6- and 7-vertices survive, 
7J5 = n-j. It is a consequence of topological constraints (1) and (2), because 
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Figure 4: Equilibrium distribution P(n) at various f3, for N = 65025. 
2 Ni = 5 ns + 6 n6 + 7 n7 

so that 

5 n 5 + 6 n 6 + 7 n7 = 6 (n 5 + n6 + n 7 ) 

hence 

n5 — n-j = 

Baby universes also disappear at low temperature, as, otherwise, they would increase curvature energy. 
However, as can be seen in simulations, there is no sign of phase transition between high and low 
temperature phases. This is also suggested by dimensional analysis, as square curvature is not relevant 
in the infrared limit. Moreover, the same conclusion arises in the matrix model mentioned above. 

3 Dynamics 

Although static properties are rather simple, dynamics of the model reveals a very rich structure. 

3.1 Finite temperature 
3.1.1 The rules 

The rules of dynamical evolution are based on the energy (3) and use so-called Alexander moves [17], 
more precisely T\ moves, which consist in flipping a link in a triangulation (figure 5). Topology and 
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Figure 5: T\ move. 



number of vertices N are conserved during evolution. It has been shown that T\ moves allow ergodic 
explorations of the space of triangulations [12]. 

One elementary Monte-Carlo step of evolution is done as follows: 

- Random choice of a link in triangulation T; 

- Flipping of the link with probability (Glauber type) 

where AE = E(7~') — E(T). The proposed triangulation 7~' is rejected if it contains vertices with 
qi < 3, self-energy diagrams or tadpoles. At zero temperature ((3 = +oo), the transition probability 
becomes 

w{T -» T') = if AE > 

= - if AE = 
2 

= 1 if AE < 

In the following, times are expressed in unit of sweeps, i.e Ni successive elementary Monte-Carlo steps. 
3.1.2 Quench to low temperature 

Sherrington et al. [4] have studied the behaviour of the system after a quench from initial infinite 
temperature (disordered) state to low - but finite - temperature. Their main result is the emergence 
of glassy behaviour. It may looks rather surprising as the model presents neither frustration nor 
quenched disorder. There are two regimes in the evolution, depending on the temperature of the 
quench. For high temperatures, equilibrium is reached after a rapid relaxation, with equilibration 
times independent of the temperature. 

On the other hand, for low temperatures, the system evolves in two steps. First, a rapid relaxation 
that corresponds to annihilation of all topological defects except 5- and 7-vertices, with fast decay of 
energy. Then, the system becomes glassy with equilibration times exponentially growing with inverse 
temperature (Arrhenius behaviour). In this second period, the energy reaches a plateau, while two 
times correlation functions of the curvature (local energy) develop shoulders and loose invariance in 
time translation. Glassy behaviour is the result of creations, annihilations and diffusions of 5- and 
7-vertices. Simplified reaction-diffusion models give similar results [18]. The authors conclude that 
the model displays properties characterizing strong glass forming systems [19]. 
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4 Zero temperature and freezing 



The particular case of deep quench to zero temperature with infinite cooling rate [20] is now described. 
This is the original contribution of this paper. 

4.1 Energy 

The evolution of the energy E(t) after a quench from infinite temperature state to zero temperature 
is plotted in figure 6 for several sizes of the system. Immediately after the quench, there is a very fast 
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Figure 6: Evolution of energy after a quench to zero temperature, for N = 2500, 6400 and 10000. 
Each curve corresponds to one run. 

decay of the energy, followed by a slowing down period where E(t) reaches a serie of plateaux. At 
the beginning, this looks like stairs. But as time goes, the plateaux become broader and broader and 
eventually "infinite". It strongly depends on the size of the system. However, simple contemplation 
of the graph of E(t) is not enough to know whether energy really stops or not. E(t) could as well 
decrease to lower plateaux, and eventually to zero, after very long times unreachable by reasonable 
computer simulations. However, as will be shown in the following, it is very likely that, after a deep 
quench, the system gets trapped in metastable states. 
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The final value reached by the energy follows a sharp distribution centered on a mean value 
Ef, according to initial disordered states. For each data point in figure 7, 16 independent runs are 
performed with 1000 heating sweeps at infinite temperature followed by sufficient sweeps at zero 
temperature (up to 5.10 5 sweeps for N = 12500) - depending on the time needed by energy to get 
frozen. The first important point is that Ef is strictly greater than zero: the system is trapped in 
metastable states and never goes to its ground state. In other words, the energy freezes, in a time that 
depends on the size of the system. This can be seen by looking at the evolution of P(n), which shows 
a very fast invasion of surfaces by 6-vertices. Both descriptions - in terms of energy and in terms of 
defects - are equivalent as energy and number of defects ndefects are simply related, if defects are 5- 
and 7-vertices only: 

E = n de f ects 

The second important point is that Ef/N is constant for large N, as can be seen in figure 7. A 
power law fit gives a limiting value lmiAr^oo Ef/N = 0.1342(4). So, mean frozen energy per site is 
a thermodynamical characteristic of the model. Alternatively, a finite fraction - ~ 13% - of defects 
remains in frozen states. 
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Figure 7: Frozen energy per site Ef/N, versus jj. The dashed line indicates a power law fit. 
4.2 T\ Moves 

At zero temperature, besides 6-vertices, only 5- and 7-vertices survive. Evolution of the system can 
be seen as dynamics of such defects through the six possible T\ moves shown in figure 8. Moves d), e) 
and /) have negative energy balance (AE < 0) and correspond to annihilation of defects. The others 
- a), b) and c) - can be viewed as diffusion of defects on the surface. In this second category, moves 
a) and b) - called flippers in the following - perform local rotations of 5-5 or 7-7 pairs. On the other 
hand, move c) is a true diffusion of a dipole - i.e a pair of 5-7. In order to understand the evolution of 
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Figure 8: Xi moves allowed at zero temperature. When there are no dipoles (5-7), only the two first 
moves (in the box) can occur. 

the system after a deep quench, it is important to know the possible ways of instantaneous evolution 
of the system at a given time. This depends on the repartition of T\ moves and gives a partial view 
on the energy landscape [21]. More precisely, at a given time t, it is possible to group the links of the 
surface into three categories: the first one (I) contains links that cannot be moved because energy 
would be increased; the second one (II) contains links whose moves make energy strictly lower if they 
are flipped; and the third one (III) contains links that let energy unchanged. Two subcategories 
of (III) can be defined: (Ilia) contains links in (III) correponding to dipole diffusion, and (I I lb) 
contains links in (III) corresponding to flippers. Proportions of links in (I), (II), (Ilia) and (I I lb) 
at time t after a deep quench are respectively written wi(t), wn(t), wuj a (t) and wnib(t)- They are 
plotted for different sizes in figure 9, obtained by averaging about 8 independent runs for each size. It 
is not surprising to see that wi(t) grows rapidly and, conversely, that wu(t) decreases rapidly: this is 
simply the evolution toward the ground state - where all links are in (I) - as expected for a system at 
zero temperature. However, after this first rapid evolution - corresponding to rapid decay of energy - 
there is a slowing down period where wji(t) is very small, meaning that the system can hardly lower 
its energy. The reason is that the system is almost unable to found local configurations that could 
lead to d), e) or /) moves. On the other hand, wju a (t) is strictly greater than zero, so that diffusions 
of dipoles still occur. Through these diffusions, dipoles can eventually get close to other defects and 
give local configurations of type (//). However, creation of such configurations by this mechanism is 
rather unlikely because dipoles are rare in this phase (see figure 10). As a consequence, the dynamics 
is very slow. This mechanism is confirmed by carefully looking at wu(t) and wju a (t) simultaneously 
(see figure 11). Small peaks of wu(t) occur immediately after small peaks of wiu a (t), meaning that 
diffusion of pairs is the only way to create local configurations leading to energy decrease. So, the role 
of (Ilia) is crucial in the slowing down period. 

However, after some time, wji and wuj a become equal to zero: this coincides with energy freezing 
and removal of dipoles. But the dynamics is not really stopped, because during the slowing down 
period described above, wnib reaches a plateau with small but finite value - independant of the size - 
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Figure 9: Probabilities of flips wj, wjj, wfjj and w\ n versus time for iV=2500, 4900 and 10000. 
Curves corresponding to different sizes are not distinguishable. 

so that flippers still occur even after removal of dipoles and energy freezing. 

At this point, the problem is to know whether flippers eventually allow diffusion of defects, through 
some collective moves, that could lead, after a long time, to creation of dipoles. In this case, energy 
could decrease again, and freezing of energy would be an illusion due to short time measurements. 

4.3 Geometric description 

In order to elucidate the role of flippers, it is necessary to have spatial informations on surfaces when 
removal of dipoles and energy freezing - at least apparently - occur. So, it is very instructive to draw 
the surfaces. As can be seen in figure 12, not only do the defects of different kind - i.e 5- and 7-vertices 
- separate, but, conversely, defects of the same kind group together into clusters. More precisely, frozen 
surfaces look like a sea of 6-vertices with isolated clusters of defects containing either 5-vertices or 
7-vertices. Figures 13 and 14 show the size distributions of these clusters, measured on 16 independent 
frozen surfaces obtained after 1000 sweeps of thermalization followed by sufficient sweeps (up to 4.10 5 ) 
at zero temperature. It should be noticed that the sizes do not exceed few vertices. Moreover, the 
distribtution itself does not depend on the size of the surfaces, meaning that the mechanism of freezing 
is local and not influenced by long range correlations. 
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Figure 10: Fraction of dipoles versus time for N = 2500, averaged on 576 independent runs. 



4.4 Flippers but no dipole creation 

The issue is now to understand the role of flippers on a surface made of 5- and 7-clusters inside a sea 
of 6-vertices. There are two possibilities: 

- 5-vertices and 7-vertices diffuse across the whole surface and, eventually, some of them are close 
enough to form dipoles and allow decreasing of energy. In this case, evolution could be very slow, but 
the system is not really frozen; 

- 5-vertices and 7-vertices are moving but clusters are confined within spatial boundaries so that no 
dipoles are created. In this case, the system is for ever trapped in a set of metastable states. 

Arguments suggesting that flippers cannot lead to dipole creation are now given. The starting 
point is a state with clusters of different kinds (5-clusters and 7-clusters), inside a sea of 6-vertices 
(about 85% of the surface). 

The 5-clusters have the following property: if there exists a convex contour C surrounding a 5- 
cluster and containing no 7-vertex, then, at zero temperature, C will never be crossed by a 5-vertex 
coming from the cluster, i.e there will never be 5-vertices outside C. Convex contour means that the 
number of links attached to each vertex of C and going outside of C is larger or equal to the number 
of links going inside C - the links of C are not taken into account. The proof runs as follows : 
Let a region 1Z, containing only 5- and 6-vertices, be bounded by a loop C, made of links connecting 
5- and/or 6-vertices. It is supposed that there are only 6-vertices outside C. Then, moves of type a) 
are performed in 1Z : 

• The boundary C cannot be broken: consider a link I G C. The only possibility for / to be moved is 
to belong to type (I I lb) (flipper). In this case, because of the convexity of the contour, one of the 
two neighbouring 5-vertices is inevitably outside C. It would contradict the starting hypothesis: so, / 
cannot be moved. 

• Consider now a set of four vertices {a\, fa, oc2, fa} involved in a move of type a), where a±, «2 are 
5-vertices and fa, fa are 6-vertices. As shown above, the link fa- fa cannot belong to C. After the T\ 
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Figure 11: Detail of wu and w a lu for one run and N=4900. 

move, cci, c*2 are 6- vertices and 0\, 02 are 5- vertices. There are two possibilities that could lead to a 
diffusion of 5- vertices outside C: 

(i) 0\ and 02 are outside C or 

(ii) 0\ or 02 are outside C. 

Of course, the starting hypothesis prevents a\ and c*2 to be outside C. So, for the case (i), C must 
contain a\ and 0:2, but neither 0\ nor /?2- It is easy to see that such a configuration implies a non 
convex boundary C. For the case (ii), C must contain a±, c*2 and /3i (resp. 02) but not /?2 (resp. 0\). 
It also implies a non convex boundary at 0\ (resp. 02)- 

To conclude, diffusion of 5-vertices outside C is not allowed. 

Consider now the 7-clusters. Unfortunately, last arguments given for 5-clusters do not work. In 
particular, a convex boundary can be deformed by a flipper. However, it is possible to have a rather 
reliable result. Consider a 7-cluster and perform - by computer simulations - T\ moves of type 6) 
only in this cluster. Then, the pattern of the cluster evolves through flippers but, as can be seen by 
drawing surfaces, it does not extend much and never reaches any 5-cluster. These simulations take 
very few computer time because 7-clusters contain only fews vertices. So, it is possible to perform 
very extensive runs and to have great confidence level in this result. 

To summarize, it is very likely that the system gets trapped and follows a constant energy walk in 
a set of states; these states are made of 5- and 7-clusters and are subject only to flippers that do not 
destroy the cluster structure - in particular, flippers cannot create dipoles. 

5 Inherent structures and energy landscape 

The quenching procedure to zero temperature used in this work maps equilibrium states to frozen 
sets of metastable states. It is in some sense a steepest-descent dynamics to inherent states as defined 
by Stillinger and Weber [22]. This approach of glassy behaviours is based on a decomposition of the 



13 



Figure 12: A frozen surface with N = 900. 5-vertices are represented by • and 7- vertices by o. All 
others vertices are 6-vertices. Arrows indicate zooms of regions in the boxes. 

configuration space into basins, called inherent structures. Each basin contains the states mapped 
through steepest-descent to the same local minimum in the energy landscape. This method provides a 
deeper understanding of glassy behaviours in terms of fast intra-basins and slow inter-basins dynamics, 
with the definition of a configurational entropy counting the number of inherent structures with a given 
energy or, alternatively, a given free energy [23]. However, this symbolic dynamics between inherent 
structures is not necessarily relevant for of all glassy systems [24] [25]. 

For the model studied in this paper, it would be very interesting to use the Stillinger and Weber 
approach. But, as the model is discrete, the zero temperature dynamics is not a deterministic mapping 
to metastable states. Moreover, these inherent structures are degenerate, as shown in the previous 
section. So, it would require slight modifications to define unambiguously concepts like inherent 
structures, configurational entropy ... The present work is a first step in this approach as it provides a 
characterization of the metastable states reached after a quench from infinite temperature equilibrium 
states. To go further, the frozen energy per site following a quench from finite temperature to zero 
temperature has been calculated for rather small surfaces (N = 400), as preliminary result. Figure 
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Figure 13: Size distribution of 5-clusters on frozen surfaces. 



15 shows the variation of Ef(T)/N with temperature. As can be seen, this quantity is more or less 
independent of T for T > T Q with T a ~ 2, and strongly decreases for T <T Q , meaning that the system 
explores deeper valleys in the energy landscape. 



6 Conclusion 

The model of random surfaces with elastic curvature energy studied in this paper shows rather trivial 
static properties. In particular, there is no phase transition driven by curvature coupling constant - 
curvature energy simply smoothes surfaces and the whole fractal structure of baby universes disappears 
at low temperature equilibrium. 

However, the system reveals complex dynamical behaviour, though there is no frustration nor 
quenched disorder. In particular, when the system is quenched from totally disordered infinite tem- 
perature state to low temperature, glassy behaviour takes place and relaxation time becomes expo- 
nentially long as temperature is lowered. 

Glassy behaviour is dramatically strengthened for deep quench at zero temperature. In this case, 
simulations and some arguments strongly suggest that the system really freezes: after a rapid decay, 
the energy stops and the ground state is never reached. Equivalently, freezing corresponds to defects 
(5- and 7-vertices) that cannot be removed by the system. This property does not depend on the 
size of the system and should be still valid at the thermodynamical limit. In fact, even if energy 
stops, there is a remaining dynamics of the system, consisting in flippers, i.e T\ moves that locally 
rotates defects of same kind. But flippers cannot make energy evolve because of the spatial repartion 
of defects: they group together into clusters made of one kind of defects . Then, flippers cannot lead 
to diffusion of defects outside some bounded regions, and dipoles, that could lower energy, are not 
created. The system is forever trapped in a set of metastable states. Similar results are obtained 
for reaction-diffusion models [18]. However, there is still an open question: why do defects of same 
kind group together into clusters? In particular, it would be very interesting to study the evolu- 
tion of the fractal structure of baby universes during the deep quench. At infinite temperature, it 
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Figure 14: Size distribution of 7-clusters on frozen surfaces. 



is known that dn = 4. In frozen states, there are no more baby universes. But preliminary results 
suggest that the fractal dimension of frozen surfaces is significantly greater than 2 - the value for 
flat surfaces. So, a guess can be made: frozen surfaces keep memory of their initial baby universes 
structure through structure of clusters. The mechanism of clusters creation would be the following : 
each initial baby universe is more or less isolated from the rest of the surface; then, just after a quench 
to zero temperature, each baby universe tries to reach its own local ground state. But in each baby 
universe, the difference between the number of 5- vertices and the number of 7-vertices is topologically 
constrained by the boundary of the baby universe. So, after a rapid multi-relaxation period where 
boundaries of baby universes are not expected to be much modified, some regions would be left with 
an excess of 5-vertices and some others with an excess of 7-vertices. This mechanism could lead to cre- 
ation of clusters. To summarize, clusters structure would be a trace of initial baby universes structure. 

It would be very interesting to use the inherent structures approach of Stillinger and Weber for 
the system studied here, and to see wether this method can provide a better understanding of the 
glassy behaviour observed at finite temperature. The present work is in some sense a first step in this 
direction. 
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Figure 15: Frozen energy per site Ef/N versus temperature, for N = 400. 

References 

[1] F. David, in Gravitation and Quantizations, (Les Houches Summer School, Session LVII, 1992) 
679-750, edited by B. Julia and J. Zinn- Justin (Elsevier, Amsterdam, 1995); 
J. Ambj0rn, in Fluctuating Geometries and Field Theory, (Les Houches Summer School, Session 
LXII, 1994) 77-193, edited by F. David, P. Ginsparg, and J. Zinn-Justin (Elsevier, Amsterdam, 
1996) ; 

P. Di Francesco, P. Ginsparg and J. Zinn-Justin, Phys. Rep. 254 (1995), 1; 

J. Ambj0rn, B. Durhuus, and T. Jonsson, Quantum Geometry, Cambridge University Press, 

(1997). 

[2] Statistical Mechanics of Membranes and Surfaces, Proc. of the Fifth Jerusalem Winter School 
for Theoretical Physics, 1987/1988, eds. D. Nelson, T. Piran, and S. Weinberg, World Scientific, 
Singapore (1989). 

[3] M.J. Bowick, S.M. Catterall and G. Thorleifsson, Phys. Lett. B391 (1997), 305. 

[4] T. Aste and D. Sherrington, J. Phys. A32 (1999), 7049; 
L. Davison and D. Sherrington, J. Phys. A33 (2000), 8615; 

L. Davison, D. Sherrington, J. P. Garrahan and A. Buhot, J. Phys. A34 (2001), 5147; 
A. Buhot, J.P Garrahan and D. Sherrington, J. Phys. A36 (2003), 307. 

[5] J. Frolich, in Lectures Notes in Physics, Vol. 216, ed. L. Garrido, Springer (1985); 
F. David, Nucl. Phys. B257 (1985), 45; 
V.A. Kazakov, Phys. Lett. B150 (1985), 182. 

[6] C. Itzykson and J-M. Drouffe, Statistical Field Theory, Vol. 2, Cambridge University Press, (1990). 

[7] T. Regge, Nuovo Cimento 19 (1961), 558. 

[8] W.T. Tutte, Can. J. Math. 14 (1962), 21. 

[9] E. Brezin, C. Itzykson, G. Parisi and J.B. Zuber, Comm. Math. Phys. 59 (1978), 35. 

[10] S. Jain and S.D. Mathur, Phys. Lett. B286 (1992), 239. 



17 



[11] H. Kawai, N. Kawamoto, T. Mogami and Y. Watabiki, Phys. Lett. B306 (1993), 19. 

[12] D.V. Boulatov, V.A. Kazakov, I. K. Kostov and A. A. Migdal, Nucl. Phys. B275 (1986), 641. 

[13] C. Godreche, I. Kostov and I. Yekutieli, Phys. Rev. Lett. 69 (1992), 2674. 

[14] D.V. Boulatov and V.A. Kazakov, Phys. Lett. B184 (1987), 247. 

[15] N. Tsuda and T. Yukawa, Phys. Lett. B305 (1993), 223; 

T. Yukawa, N. Tsuda and A. T. Sornborger, Nucl. Phys. B(Proc. Suppl.)30 (1993), 791; 

S. Ichinose, N. Tsuda and T. Yukawa, Nucl. Phys. B445 (1995), 295; 

S. Ichinose, N. Tsuda and T. Yukawa, Int. J. Mod. Phys. A12 (1997), 757. 

[16] V.A. Kazakov, M. Staudacher and T. Wynter, Nucl.Phys. B471 (1996), 309. 

[17] J.W. Alexander, Ann. Math. 31 (1930), 292; 

M. Gross and S. Varsted, Nucl. Phys. B378 (1992), 367. 

[18] J-M. Drouffe, private communication. 

[19] W. Kob, J. Phys. Condens. Matter 11 (1999), R85. 

[20] S. Sastry, P. G. Debenedetti and F. H. Stillinger, Nature 393 (1998), 554. 

[21] K. Broderix, K. K. Bhattacharya, A. Cavagna, A. Zippelius and I. Giardina, Phys. Rev. Lett. 85 
(2000), 5360. 

[22] F.H. Stillinger and T.A. Weber, Phys. Rev. A25 (1982), 978. 

[23] F. Ritort and P. Sollich, Advances in Physics 52 (2003), 219; 

B. Coluzzi, A. Crisanti, E. Marinari, F. Ritort and A. Rocco, Eur. Phys. J. B32 (2003), 495; 
A. Crisanti and F. Ritort, J. Phys.: Condens. Matter 14 (2002), 1381. 

[24] A. Crisanti, F. Ritort, A. Rocco and M. Sellitto, J. Chem. Phys. 113 (2000), 10615; 

A. Crisanti, F. Ritort, A. Rocco and M. Sellitto, J. Phys.: Condens. Matter 14 (2002), 1523. 

[25] G. Biroli and J. Kurchan, Phys. Rev. E64 (2001), 016101. 



18 



